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Abstract 

We present a high statistics lattice calculation of the B-meson binding energy A 
of the heavy-quark inside the pseudoscalar B-meson. Our numerical results have been 
obtained from several independent numerical simulations at /3 = 6.0, 6.2 and 6.4, and 
using, for the meson correlators, the results obtained by the APE group at the same 
values of (3. Our best estimate, obtained by combining results at different values of 
/3, is A = 180]':2o MeV. For the MS running mass, we obtain rnb{mh) = 4.15 ± 0.05 ± 
0.20 GeV, in reasonable agreement with previous determinations. The systematic error 
is the truncation of the perturbative series in the matching condition of the relevant 
operator of the Heavy Quark Effective Theory. 
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1 Introduction 



In a previous paper we have presented the first results for the binding and kinetic 
energies of the 6-quark in a i?-meson, obtained by using the lattice version of the Heavy 
Quark Effective Theory (HQET) j^. These are the non-perturbative parameters A 
and Ai, which control the ©(l/m^) and 0(l/m^) corrections to the spectroscopy and 
inclusive decays of S-mesons, and which also appear in the theoretical predictions for 
several exclusive processes [|. mb, tub and m_B* denote the masses of the 6-quark and 
of the B and B* mesons respectively. In this letter we report on a detailed study of A, 
which improves on that in ref. in two ways: 

i) We have collected a much larger set of gauge field configurations than was used in 
ref. Q , with a consequent reduction of the statistical error in the determination 
of A. Moreover, the more accurate knowledge of the heavy quark propagator has 
allowed for an improved study of its infrared behaviour (i.e. its behaviour at large 
time separations). This will be discussed in detail in the following. 

ii) We also present new results at a value of the lattice spacing which is smaller than 
those considered in ref. corresponding to (3 = 6.4. Combining the results 
at three different values of /3 (/3 = 6.0, 6.2 and 6.4) it is possible to study the 
dependence of A on the lattice spacing. Within the precision of our computations, 
we find that the results are independent of the lattice spacing. 

The presence of renormalon singularities in the pole mass [P, ^ implies that the 
most intuitive definition of the parameter A, i.e. A = uib — m^°^'^, where denotes 
the pole mass of the 6-quark, does not correspond to a physical quantity. Thus if A is to 
be introduced into phenomenological applications (which is not necessary, see below), 
alternative definitions have to be given. For example, it is possible to define A from 
the experimentally measured value of some physical quantity, for which the theoretical 
prediction depends on A. Since A does not contribute directly to the leading behaviour 
of any physical quantity, always appearing as an 0{l/mi,) correction, in order to extract 
A, it is first necessary to subtract the perturbation series for the leading term (for a 
detailed explanation see ref. Q). This series has a renormalon ambiguity, implying 
that its sum is not uniquely defined, but the subtraction is performed up to some finite 
order of perturbation theory. The physical quantity being used, and the order of the 
perturbation series which has been subtracted, constitute the definition of A 0. 

Lattice simulations provide the opportunity to compute A, defined in some suitable 
way, non-perturbatively. In ref. a definition of A was proposed which is free of renor- 
malon ambiguities, and from which the linear ultraviolet divergence (i.e. the divergence 
of 0{a~^), where a is the lattice spacing) has been subtracted non-perturbatively. The 

^In addition there is the parameter A2, which is the matrix element of the chromomagnetic operator. This 
parameter gives the leading term in the hyperfine splitting (m^. — m^), and so can be estimated directly 
from the experimental data. 

^ Notice that such definitions of A retain some dependence on the mass of the heavy quark. 
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necessity for non-perturbative subtractions of power divergences has been stressed in 
ref. 0. This definition of A, which will be explained in section |2|, is obtained after 
imposing a non-perturbative renormalisation condition on the heavy quark propaga- 
tor in the Landau gauge. Even though, in practice, the value of A defined in this 
way is obtained using lattice simulations, it is in fact independent of the method of 
regularisation. 

Although it may be convenient to try to introduce a parameter A, which is of 
O(Aqcd) and is free of renormalon ambiguities, it is not necessary to do so. One 
can either relate two or more physical quantities directly to the required precision ^, 
or determine physical quantities from non-perturbative computations, such as lattice 
simulations Below, as an example of the computation of a physical quantity, we 
present the results for the 6-quark mass (nih) renormalized in the MS scheme at the 
scale ji = rrif). The main uncertainty in this quantity comes from the truncation at 
0(as) of the perturbative factor necessary to match the quark mass computed in the 
HQET to mft. We briefly discuss how to reduce this error which we currently estimate 
to be about 200 Me^^g. 

Any deflnition of A which is of O(Aqcd) and is independent of any renormalization 
scale, necessarily involves Green functions deflned at large distances. Although we 
have introduced such a definition in ref. |^], which is free of renormalon singularities, 
we cannot prove that non-perturbative confining effects do not invalidate the proposed 
procedure. An important aspect of this paper is to study the infrared behaviour of the 
heavy quark propagator in order to search for such effects. They are not visible within 
the precision of our lattice calculations, up to distances of about 1.5 fm, which is the 
largest distance which we can reach. 

The plan of the paper is the following. In sec. ^ we recall the relevant formulae for 
our definition of A Q; in sec. ^ we describe the numerical calculation of A and nib and 
discuss the main results of this study; in the conclusions we present the outlook for 
future developments and applications of the method discussed in this paper. 



2 Non-perturbative definition of A 

In this section we define the non-perturbative renormalisation prescription which we 
will use to calculate the "physical" value of A H]. A detailed discussion of our procedure 
has been presented in our previous papers ||l[, |^, and so we only exhibit here those 
formulae which are necessary to the understanding of the numerical results. 
Consider the following Lagrangian of the lattice formulation of the HQET: 

CcS = ? — (h{x)D4h{x) + 6mh{x)h{x)), (1) 

1 + a dm V ^ 

•^This is done implicitly when A is defined using one physical quantity, and used in the prediction for a 
second one. 
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where the factor 1/(1 + a6m) has been introduced to ensure the correct normalization 
of the heavy quark field h, and the covariant derivative is defined by 



f{t) = l/a (/(t)-[/t(t-a)/(t 



(2) 



Lattice computations are frequently performed by choosing the bare-mass {6m) to be 
zero. We propose instead to choose a value of the bare mass by imposing a renormalisa- 
tion condition on the heavy quark propagator. With our prescription, the heavy-quark 
propagator can be made finite, up to the standard logarithmic divergences associated 
with the renormalization of the heavy-quark field. 

The mass counter-term can be chosen in many different ways. One possibility is to 
define it by studying the behaviour of heavy-quark propagator S{x, t) in some gauge 
(we will use the Landau gauge) at large values of the time: 



— 5m = 



ln(l -|- a 6m) 



lim — In 



Tr 



{S{x, 



t + a] 



Tr 5(f, t) 



(3) 



where the traces are over the colour quantum numbers, and we have assumed that the 
large time limit of the ratio in eq. (^) exists. Our numerical results, to be discussed in 
detail in the next section, support the validity of this assumption and the use of eq. (^ 
is our preferred definition of 6rn. In words, (^) is the condition that the subtracted 
heavy quark propagator has no exponential fall (or growth) at large times. If, instead, 
the bare mass 6m is chosen to be zero, loop corrections generate a mass of 0(1 /a), 
which leads to an exponential behaviour in time. 
A possible alternative definition is given by 



- 6m{t* 



In 



Tr 



(5(f , t* 



+ a, 



TV S(x, t* 



(4) 



for some chosen time t*. The corresponding definition of A, see eq. (0) below, will 
clearly depend on the choice of t*: t* parametrizes the renormalisation prescription 
dependence and can be considered as the renormalisation point in coordinate space. 
For physical matrix elements the dependence on t* is eliminated in the matching of the 
QCD action and operators with those of the HQET. The use of the propagator at small 
times, t*AQCD ^ 1, to define 6m{t*), and hence A(i*), does not require any assumption 
about the large time behaviour of the heavy quark propagator and in section ^ we will 
also present the results for A(t*) obtained in this way. However A(t*) contains terms 
of 0{l/t*), and is therefore not a parameter of O(Aqcd)- 

The divergent parts of the counterterms, either 6rn or 6m{t*), are gauge invariant, 
in spite of the fact that they are calculated from the heavy quark propagator in a fixed 
gauge. The argument goes as follows. The linear divergence is eliminated from any 
correlation function, i.e. for any external state, by subtracting from the action a term 
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proportional to the gauge-invariant operator hh. Since in this way one eliminates all 
divergences both for quark and hadron external states, the coefficient of the mixing 
has to be gauge- invariant. This must be true also for the finite non-perturbative con- 
tributions which arise from the coefficient of the linear divergence. These contributions 
are generated by terms which are exponentially small in the inverse coupling constant, 
see ref. |^, 0. As for the finite parts, the gauge dependence, when it is present, is 
eliminated from physical quantities up to the order in perturbation theory at which 
the calculation is performed, see e.g. the case of the running mass discussed below. 

In order to obtain the renormalized A, in addition to the mass counter-term 5rn, we 
need to compute, the "bare" binding energy £. In the lattice HQET £, which diverges 
linearly in 1/a, is computed from the two-point correlation function 

C{t) = J2{0\J{x,t)j\0,0)\0), (5) 

X 

where J (J^), is any interpolating operator which can annihilate (create) a S-meson ^. 
For a sufficiently large Euclidean time t, 

C{t) exp(-a) (6) 

where Z is a constant. 

We now define the renormalised binding energy by 

A = £ - dm, (7) 

and a "subtracted" pole mass by 

rub = TUB — £ + 5m . (8) 

A and mf, defined in this way are free of both renormalon ambiguities and linear 
divergences. Similar formulae hold in the case when 5m{t*) is used instead to define 
A{t*) and mb{t*). 



3 Numerical implementation of the renormali- 
sation procedure 

As explained in the previous section, the determination of A requires the computation 
of the heavy quark propagator in a fixed gauge (in order to obtain the subtracted 
operator), as well as the evaluation of the two-point correlation function (^). We 
have obtained our results using three independent numerical simulations, whose main 
parameters are given in table |l|. 

^In this paper we only consider pseudoscalar B-mesons, but an analogous discussion holds for other beauty 
hadrons, such as the A^. 



4 



simulation 


volume 


/3 


Number of configurations 


set A 


24^ X 40 


6.0 


500 


set B 


24^ X 40 


6.2 


300 


set C 


24^ X 40 


6.4 


300 



Table 1: Parameters of the numerical simulations, the results of which have been used for 
the present study. 



Our best value of the subtracted binding energy, K = 8 — 5rn, has been determined 
by combining the values of 5rn obtained using set A . set B and set C with the calcu- 
lation of £ performed by the APE collaboration at /3 = 6.0, 6.2 [0] and 6.4 which will 



be published elsewhere yj. £ was determined using the SW-Clover fermion action |1C] 
for the light quarks in the quenched approximation. The calculations were performed 
at several values of the masses of the light quarks, so that extrapolations to the chi- 
ral limit are possible. 5rn {5rn{t*)) has been determined by using either the standard 
definition of the lattice heavy-quark propagator or its improved (up to and including 
0{a) terms) version, the explicit espression of which can be found in ref. [||]. Only the 
results obtained by using the improved heavy-quark propagator are reported in the 
following. 

All the errors have been computed with the jacknife method by decimating 10 
configurations at a time for set A , and 6 configurations at a time for set B and set C . 
The error on £ was computed with the jacknife method also and we refer the reader 
to refs. [^, ^] for details. 

3.1 Determination of the residual mass drn 

To determine the residual mass, we have to compute the effective mass 

of the heavy quark propagator. Suit), defined as 

Snit) = ^J2{TT[S{t,x)]), (10) 

X 

where the trace is over the colour indices and V denotes the spatial volume of our 
lattice. It is this averaged propagator Snit), which has been used in the computations 
below; the averaging over all spatial points at time t reduces the statistical errors 
enormously. 

In fig. 1^ we present the values of 6m{t) for the improved propagator as a function of 
t/a. The figure shows the very high accuracy which can be reached in the determination 
of the effective mass, i.e. in the determination of the logarithmic derivative (p[) of the 
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Figure 1: Effective mass of the improved heavy-quark propagator Snit), obtained at (3 — 6.0 
using 500 configurations, and at (5 = 6.2 and 6.4 using 300 configurations, as a function 
of the time. The curves represent fits of the numerical results to the expression given in 
eq. (12). 
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heavy quark propagator. This will allow us to make a careful study of the infrared 
behaviour of the quark propagator. 

As shown in ref. [|^, the effective mass is indistinguishable in the improved and 
unimproved cases, for t/a> 4-5. Thus, in order to minimize lattice artefacts, we have 
only used the results obtained for t/a > 6. 

Inspired by the results of one-loop perturbation theory Q, we made a fit to 5rn{t) 
using the expression 

a 5rn{t) = a 5rn + "f — . (11) 

In order to mimic higher order effects, we have also used different expressions to fit 
6fn{t), e.g. 

a5m(t) = a 5m + y In (^^) (12) 



or 



and changed the interval of the fits in order to check the stability of the determina- 
tion of 6fn. Finally, in order to monitor possible confining effects in the heavy-quark 
propagator, we have also tried a fit of the form 

a Srn{t) = a 5rn + J ln(^^^^^ + pt , (14) 

where the last term in eq. (|14|) will be discussed later on. In eqs. (pTl)-(p!^), 5m, 
7, . . . , 7", C and p are free parameters of the fit. 

From the different results obtained by varying the fitting functions (|ri|)-(p!3|) and 
the time intervals, we quote 

adfn = 0.526 ± 0.003 ± 0.006 at /? = 6.0 (15) 

a6fn = 0.458 ± 0.005 ± 0.004 at /3 = 6.2 (16) 

adrn = 0.407 ± 0.006 ± 0.006 at /? = 6.4 (17) 

where in all cases the first error is statistical, and the second is an estimate of the 
systematic uncertainty, based on the spread of results obtained using different time 
intervals and fitting functions. Notice the improvement of the above figures, when 
compared to our old results a 6m = 0.521 it 0.006 ± 0.010 at /? = 6.0 and a 6m = 
0.445 ± 0.008 ± 0.010 at (3 = 6.2. The improvement is due to the increased statistics, 
which also allowed us to study the heavy quark propagator at larger time distances, 
thus reducing also the systematic error. 
In order to obtain A we have used: 

• 6fn from eqs. ([T5|)-(p!7|); 

• the SW-Clover determination of £ of the APE collaboration, a£ = 0.610 ± 0.010 
atp = 6.0, a£ = 0.520±0.010 at /3 = 6.2 i] and a£ = 0.460±0.007 at /3 = 6.4 [|; 
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• a'^ip = 6.0) = 2.0 ±0.2 GeV, a-i(/? = 6.2) = 2.9 ±0.3 GeV and a'^iP = 6.4) = 
3.8 ± 0.3 GeV. A significant contribution to the final error in A comes from the 
calibration of the lattice spacing in quenched simulations which typically has an 
uncertainty of 0(10%), depending on the physical quantity which is used to set 
the scale. We take the results given above as a fair representation of the spread 
of values of a~^. 

We then find 

A = £■ - (5m = 170 ± 30 MeV at /? = 6.0 (18) 

A = £■ - (5m = 180 ± 40 MeV at /? = 6.2 (19) 

A = £: - (5m = 200 ± 40 MeV at /? = 6.4 (20) 

where the statistical errors have been combined in quadrature with those due to the 
uncertainty in the lattice spacing. 



Within the uncertainties, the results in eqs. (18)-(pO|) are compatible with the 
expected independence of A of the lattice spacing. Given the intrinsic uncertainty in 
the value of the lattice spacing in quenched simulations it is difficult however, to check 
this more precisely. Certainly our results exclude large 0{a) effects. Assuming that A 
is indeed constant, we combined together the results of eqs. (|T8|)"(pO[) to obtain 



A = (ISOtll) MeV , (21) 

where in the final error we have included also the discretization error which has been 
conservatively estimated to be of the order of +20 MeV in ref. [|| . The estimate of the 
discretization error was obtained by comparing the value of £ obtained by using the 
standard Wilson action and the SW-Clover action at the same values of (3. 

It would be interesting to compare this value of A with other determinations, ob- 
tained by using different methods, e.g. the recent result obtained from the semileptonic 
width in ref. ||ll|. The relation between our definition and other ones, provided they 
are implemented consistently and are free of renormalon ambiguities, can be formally 
derived in perturbation theory. Even in the case where the perturbative series which 
connects two different definitions of A is precisely determined, this would not reduce 
the error in the prediction for physical quantities, such as nih- Indeed, as discussed 
below, for m;, the error of about 200 MeV is dominated by the perturbative series 
which matches QCD to the effective theory, and not by the error in A. 

The determination of the mass counter-term 6rn(t*) at fixed t = t* , requires no 
fitting, and the results obtained using set A-set C are presented in table ^. Only the 
results obtained in a range of t* corresponding to values of n/t* between 1.0 and 2.5 
GeV are reported in the table. The results for A(t*) obtained using the same values of £ 
and of a^^ as before, are also reported in the same table. They show that the values of 
the binding energy at the same value of the renormalisation scale tt /t* (but at different 
values of the lattice spacing) are very consistent, implying that lattice artefacts are 
rather small. 
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(3 


t*/a 


aSm{t*) 


Tv/t* (GeV) 


A(r) (MeV) 


6.0 


3 


0.3694(7) 


2.1 ±0.2 


480 ± 20 ± 50 




4 


0.4004(8) 


1.6 ±0.2 


420 ± 20 ± 40 




5 


0.4207(8) 


1.3 ±0.1 


380 ± 20 ± 40 




6 


0.4356(8) 


1.0 ±0.1 


350 ± 20 ± 40 


6.2 


4 


0.3478(6) 


2.3 ±0.2 


500 ± 20 ± 50 




5 


0.3643(6) 


1.8 ±0.2 


450 ± 20 ± 50 




6 


0.3765(10) 


1.5 ±0.2 


420 ± 20 ± 40 




7 


0.3868(12) 


1.3 ±0.1 


390 ± 20 ± 40 




8 


0.3949(15) 


1.1 ±0.1 


360 ± 20 ± 40 


6.4 


5 


0.3252(4) 


2.4 ±0.2 


510 ± 30 ± 40 




6 


0.3355(4) 


2.0 ±0.2 


470 ± 30 ± 40 




7 


0.3435(5) 


1.7 ±0.1 


440 ± 30 ± 30 




8 


0.3499(7) 


1.5 ±0.1 


420 ± 30 ± 30 




9 


0.3561(7) 


1.3 ±0.1 


390 ± 30 ± 30 




10 


0.3604(13) 


1.2 ±0.1 


380 ± 30 ± 30 




11 


0.3669(18) 


1.1 ±0.1 


350 ± 30 ± 30 



Tabic 2: Results for A(t*) = S — 5m{t*) for different normalization times t* , using the 
results from set A-set C . The first error on A(t*) is obtained by combining the errors on £ 
and 5m{t*) in quadrature; the second error (and the error on n/t*) comes from the calibration 
of the lattice spacing. 
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We conclude this subsection with a brief discussion of the infrared behaviour of the 
heavy quark propagator. Due to confinement effects, it is not clear that the logarithmic 
derivative (^) has a finite limit as t ^ oo. Indeed, if ^ goes to a constant at large 
time distances, then the propagator in the effective theory behaves like an on-shell 
propagator of a particle of mass 5rn, Snit) ~ exp(—6rnt). This contradicts the naive 
expectation that the Fourier transform of the propagator of a confined object should not 
have poles, as was found, for example, in lattice computations of the gluon propagator 
[12|. In the case of the heavy-quark propagator, it is also possible that the "effective" 
mass 5rn{t) does not tend to a constant at large time-distances, and that the propagator 
decreases faster than an exponential. In order to monitor this, we have allowed for a 
linear term in t in the effective mass and have fitted (P) to (|l4|). We found that the 
slope p is compatible with zero |^ for any fitting interval, and at all the values of j3 
considered in this study. Morevover the value 5rri obtained by fitting 5m{t) to eq. (12) 
or to eq. ([T^ ) are equal within the errors. This implies that no sizeable variation of 
the effective mass, due to confinement effects, is observed down to rather low scales 
(as small as vr/lGa ~ 390 MeV, or distances up to 1.6 fm, at /? = 6.0), which belong to 
the non-perturbative region. 



3.2 The MS mass of the 6-quark 

Although it is convenient and conventional to introduce parameters such as A when 
studying the 0{l/inh) corrections to quantities in heavy quark physics, it is not nec- 
essary. One can instead derive physical quantities from results obtained using lattice 
simulations. In this subsection we present the results for the MS mass of the 6-quark, 
from simulations in the Heavy Quark Effective Theory The evaluation of a short 
distance mass does not require any assumptions about the long-distance behaviour 
of the quark propagator, indeed it can be obtained from computations of the meson 
propagator only. In ref. it was shown that, at one-loop order, 

m,= {Ms-£ + aM)^){l-'-^), (22) 

where X is the coefficient of the one-loop mass term, obtained in the lattice HQET 
with zero bare mass. The second factor on the right-hand side of eq. (^) is that which 
relates the pole mass to rrib. The perturbative series of linearly divergent terms (of 
which only the one- loop term is explicitly exhibited in eq. (^)) removes the linear 
divergence form £, but has a renormalon ambiguity of O(Aqcd)- This ambiguity is 
cancelled by that in the second factor, i.e the series relating nib to the pole mass. Thus 
nib h^-s neither a linear divergence not a renormalon ambiguity. In evaluating nib we 
take the experimental value of Mb (5.278 GeV) and the value of £ as measured in the 
lattice simulations of the APE collaboration. as{a) is taken in the range 0.13 and 0.18. 

^ For example p = 0.002 ± 0.002 at /3 = 6.0 by fitting in the interval 6<t/a< 16. 
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We also evaluate mi, using 



nib = mb{t*) X Cmit*) , 



where 



and 



mbit*) = Mb - Ait*) 
4as(mfe) 



Mb - £ + 5m{t*) , 



Cmit*) = 1 



Svr 



^ ( Smit*) 
rrib \ 



a. 



a J 



(23) 



(24) 



(25) 



Eq. ( psD is equivalent to (|22| ) up to one-loop order and neglecting terms of 0(l/mb) 
in the quark mass. Numerical differences in the results obtained in these two ways are 
a partial measure of the systematic uncertainties. The values of admit*) used in the 
analysis are taken from eqs. (|^) — (^) and table |2|. 

To obtain a distribution of results, we varied f, Aqcd and admit*) according to 
a gaussian distribution; a~^ was varied with flat distribution within its error, while 
as (a) was written in terms of the leading quenched expression of the running coupling 
constant, evaluated at the scale vr/a, with Aqcd distributed according to a flat distri- 
bution of width a and such that = 0.13 for Aqcd — cr and as = 0.18 for Aqcd + c^- 
The resulting distribution is a pseudo-gaussian [ffl] . From the width of the distribution 



we estimate the average value and error on mb. The results obtained using eqs. (23) 
and (22) are given in table ^. 



MS running 6-quark mass, mbimb) 


a 


Eq. (24) t* = 3 - 6 


Eq. (24) t* = oo 


Eq. (22) 


nibimb) (GeV) 


6.0 


4.20(7) 


4.18(7) 


4.22(7) 


4.20(7) 


6.2 


4.13(9) 


4.10(9) 


4.15(8) 


4.13(9) 


6.4 


4.07(9) 


4.04(9) 


4.09(8) 


4.07(9) 


AVERAGE 


4.15 ± 0.05 ±0.05 



Table 3: The MS running mass of the b quark. The first error is statistical and the second 
systematic. 

From the results in the table, we estimate frib = 4.15 it 0.05 ± 0.05 GeV, where 
the second error is the systematic error, evaluated by comparing the results obtained 
using the different formulae given in eqs. (|2^) and (p^. This error, which is of about 
50 MeV, can be interpreted as being due to higher order corrections in a^. It can be 
compared with the estimate of the effects of higher order corrections, derived from the 
the bubble-resummation approximation, which is about 200 MeV Q. We believe that 
an estimate of 200 MeV is a more realistic one | and for this reason we quote as our 
final result 



frib = 4.15 ± 0.05 ± 0.20 GeV . 



(26) 



'The ignorance of higher order corrections also implies an error of this order in any approach. 
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4 Conclusions 



In this study we have presented a high statistics calculation of the binding energy of 
the 6-quark in a i?-meson. The improved statistics, with respect to ref. |l| has allowed 
a better control of the infrared behaviour of the heavy quark propagator. There is no 
point in further increasing the statistics in the calculation of 6rn, since the precision on 
A is mainly limited by the calibration of the lattice spacing in quenched simulations 
and by the precision in the determination of the bare binding energy £. It is difficult 
to imagine a reduction of the uncertainty in the calibration of the lattice spacing in 
the quenched case, while the error on the bare binding energy is likely to be reduced 
in the near future. 

The main error in the evaluation of the MS 6-quark mass (~ 200 MeV) comes 
from the unknown higher order terms in the perturbative matching of the HQET to 
the full theory j^. This error can be reduced by computing perturbative corrections 
at two- loop order. This requires the calculation of the heavy quark propagator in the 
continuum and of the propagator in the lattice effective theory at O(a^). The contin- 
uum calculation has already been performed [13|. We believe that the best approach 
to the two-loop (and higher order) calculations in the lattice effective theory is the one 
based on the Langevin equation being developed by the authors of ref. Based 
on the summation of light-quark bubble graphs [l^], we estimate that knowledge 
of the matching factor at two-loops would reduce the theoretical uncertainty to about 
100 MeV. 
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